Magnetic Reconnection by a Self-Retreating X-Line 



OO 
O 

o 

o 
O 



Oh 
I 

O 



(N 
> 

00 

o 

OO 

o 



X 



M. Oka,i'B M. Fujimoto,^ T. K. M. Nakamura,^ I. Shinohara,^ and K.-I. Nishikawa^ 

^ Center for Space Plasma and Aeronomic Research, 
University of Alabama in Huntsville, Alabama, 35805, USA 
'"^ Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, Sagamihara, Japan 

(Dated: October 24, 2008) 

Particle-in-cell (PIC) simulations of collisionless magnetic reconnection are performed to study 
asymmetric reconnection in which an outflow is blocked by a hard wall while leaving sufficiently 
large room for the outflow of the opposite direction. This condition leads to a slow, roughly con- 
stant motion of the diffusion region away from the wall, the so-called 'X-line retreat'. The typical 
retreat speed is ~0.1 times the Alfven speed. At the diffusion region, ion flow pattern shows strong 
asymmetry and the ion stagnation point and the X-line are not collocated. A surprise, however, is 
that the reconnection rate remains the same unaffected by the retreat motion. 



PACS numbers: 52.35.Vd, 52.65.Rr, 94.30.cp, 96.60.qe 

Magnetic reconnection triggers many explosive phe- 
nomena in laboratory and astrophysical plasma. Impor- 
tance lies in the diffusion region where the MHD breaks 
down and the kinetic scale comes in. It is a scientific chal- 
lenge to understand the structure of the diffusion region 
and its role in controlling the reconnection rate. Owing 
to the growing power of computer capabilities, full parti- 
cle simulations has become popular these days. It is now 
clear that the out-of-plane current remains localized and 
the reconnection rate remains fast while electrons form a 
high- velocity jet that extends large distances downstream 
from the X-line [3, [3, ^, 41 • The region with the local- 
ized current is termed 'inner diffusion region' whereas the 
rest of the diffusion region including the elongated elec- 
tron current layer is termed 'outer diffusion region'. In 
general, simulations generate magnetic reconnection that 
is symmetric in both inflow and outflow directions. 

Actual magnetic reconnection, however, occurs in a 
more complicated situation. During a magnetic recon- 
nection in the Earth's magneto-tail, an earthward out- 
flow soon collides with the dipole field while the tailward 
outflow is directed to the interplanetary space. As for 
the solar flares that take place above the area of emerg- 
ing fluxrope, the upward outflow eventually merge into 
the solar wind while the downward outflow collides with 
the magnetic obstacle in the lower part of the corona. In 
these cases, magnetic reconnection may not be symmet- 
ric about the X-line. 

In this letter, we investigate the consequences of block- 
ing one side of the outflow. A hard wall is set up just 
ahead of an outflow while leaving a sufficiently large 
room for the outflow of the opposite direction. As a re- 
sult, reconnected magnetic fields are piled-up against the 
wall. By the time the pile-up region reaches the diffu- 
sion region, the X-line starts to move away from the wall 
('X-line retreat'). Under such a strong influence of the 
boundary condition, the structure of the diffusion region 
is modulated. In the past, the X-line motion was consid- 
ered by an analytical treatment fsj , MHD simulations Q , 



and a full particle simulation 'T] , although these work did 
not shed light on the modulated structure of the diffusion 
region. This paper will be the first to report the struc- 
ture of the diffusion region in a self-consistent simulation 
of asymmetric outfiow reconnection. 

We utilized the two dimensional, particle-in-cell code 
Q. The initial condition consists of a Harris cur- 
rent sheet. The anti-parallel magnetic field and the 
density are given by By = i?otanh((a; — Lx/2)/D) and 
Ncs = Nq/ cosh^{{x — Lx/2)/D), respectively, where 
Bq is the magnetic field at the inflow boundary, D is 
the half-thickness of the current sheet, Nq is the den- 
sity at the current sheet center and and Ly are the 
domain size in x and y direction, respectively. D is 
chosen to be 1.5Ai where Xi is the ion inertial length. 
The inflow, background plasma is represented by Nb = 
Nbo{1-1/ cosh^{{x-L^/2)/D)) where TVso = . 27Vo ■ The 
pressure imbalance by the non-uniform density is not im- 
portant to the results presented below. The electron to 
ion temperature ratio is set to be Te/Ti=l/5 for both the 
current sheet and the background. The frequency ratio 
cijpe/r2ce=1.5 where and f2ce are the electron plasma 
frequency and the electron cyclotron frequency, respec- 
tively. The conducting walls are used at x—0 and x—L^ 
and symmetric boundary conditions were used at y=0 
and y—Ly. At the symmetric boundary, particles are 
specularly reflected whereas fields are given by dBn/dn 
= 0, Bt = 0, d'Et/dn — 0, and En — 0, where E and 
B are the electric field and the magnetic field vectors 
and the subscripts n and t denote the normal and tan- 
gential components, respectively. Reconnection is ini- 
tiated at the distance H from the bottom (y—O) wall 
with a small magnetic island given by the vector potcn- 
tial A, = -AocM-{{x - L,/2)2 + {y - Hf}/{2X^)] 
where = Q.lBoXi. Ly is chosen to be large enough 
to simulate effectively the free boundary condition of the 
reality. 

Five simulation runs are performed to study the depen- 
dence on the initial distance H. The values of H along 
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TABLE L Initial parameters and results of each run. Lengths 
are normalized by the ion inertial length Xi, times are by the 
inverse ion cyclotron frequency and speeds are by the ion 



Alfven speed va 
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See text for details. 
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with the other initial values are compiled in Table [J The 
initial X-line is at the center of the domain for Run 4, 
and thus, it is a non-retreating case. Run 5 is an addi- 
tional case intended to study the dependence on the ion 
to electron mass ratio fi = rrii/me. In all runs, we used 
average of 64 particles in each grid cell. 276 particles per 
cell represents the unit density. 
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FIG. 1: Overview of the retreating X-line for Run 1: (a-c) 
normal electron current Jez- (d) positions of the X-line and 
the stagnation points, (e) the X-line retreat speed (black 
curve), upward and downward maximum flow speeds for ions 
(blue solid curves) and electrons (red solid curves), and the 
flow speeds at the X-line for ions (blue dashed curve) and 
electrons (red dashed curve). A stagnation point is defined 
by the flow reversal point. 



Upper panels of Figure [T] show snapshots of electron 
out-of-plane current for Run 1 . It is evident that as mag- 
netic reconnection proceeds, a downward flow is blocked 
by the bottom wall and creates the pile-up region which 
expands in time and eventually pushes the diffusion re- 
gion upward. The upward flow, on the other hand, con- 
tinues to grow in length until its leading edge reaches 
the upper wall. The choice of Ly is not sensitive to the 
results presented in this paper and thus any effect from 
the upper wall is not important. The diffusion region is 
also elongated and the final length of the electron current 
sheet becomes as long as ^^50 A.;. 

The black curve in Figure [TJi shows the X-line posi- 
tion defined by the reversal point, Bx=0. A distinct 
retreat process starts at rtdt ~85. The X-line moves 
away from the wall monotonically. Although a jump ap- 
pears at ilcit ~145 due to a magnetic island generated 
by a secondary tearing instability, the X-line motion con- 
tinues even after the upward outflow reaches the upper 
wall. The approximate time of the outflow collision with 
the wall is fldt ^180. 

Figure[T^ shows various flow speeds. While the upward 
ion flow becomes super- Alfvenic and continues to increase 
until the time of collision with the wall, the downward ion 
flow reaches ^0.3va at fldt ~80 when it is blocked by the 
bottom wall (marked by the blue arrow). The downward 
flow speed soon decreases and is followed by the upward 
motion of the X-line. In contrast, the retreat speed and 
the ion flow speed at the X-line are roughly constant at 
the same value ^O.Iva, although disturbances appear at 
^cit ^145 due to a secondary tearing. 

A striking feature can be found around the diffusion 
region. Figure [5] shows an enlarged view of the diffusion 
region for Run 1 at fldt ~115 when the upward electron 
speed is saturated. In this figure, data are accumulated 
over one gyro period. Coordinates are shifted so that 
the X-line comes to the origin. This simulation, along 
with others with different H, reveals that the ion and 
electron stagnation points are not collocated with the 
X-line. As the pile-up region expands, the lower half 
of the outer diffusion region is modulated and ion flows 
are deflected. As a result, the ion stagnation point is 
shifted downward about the distance of Xi from the X-line 
(marked by the blue horizontal arrow). This separation 
between the stagnation point and the X-line is almost 
constant in time as shown in Figure [T]i. The shift of the 
ion stagnation point leads to an upward ion flow at the 
X-line, resulting in a slow, rising motion of the X-line. 
While there is a high speed ion flow in the upper half of 
the diffusion region (the curve labeled T' in the bottom 
panel) , the typical value of the upward ion flow at the X- 
line was ~0.1va ('H'). The deflection can be recognized 
as the double peak feature in the cut of the downward 
outflow (TIF). 

In contrast to ions, electrons do not show strong deflec- 
tion and exhibit well developed outflow jets both upward 
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FIG. 2: Enlarged view of the retreating X-line for Run 1 obtained at f2cit=115. The center panel shows the electron flow in 
y direction Vey (color code) as well as ion flow directions in the x — y plane (vectors). The left panel shows quantities B^, 
Vey/y/JI, and Viy along x—Q whereas the right panel shows the electric field terms of the generalized Ohm's law along x=0. 
The bottom panel shows Viy cut along the horizontal lines in the color code. The purple and yellow arrows in the right panel 
show the outer and inner diffusion regions (DR), respectively. 



and downward, although the downward jet is limited in 
length. The typical length of the downward electron jet 
is ~5Ai and remains constant throughout the run. The 
maximum value of the electron flow speed, however, is 
quite fast and reaches 0.65vAe, as can be seen in the left 
panel. This is almost the same as the maximum speed of 
the upward electron jet. We have found that all retreat- 
ing cases show nearly the same maximum speed of the 
downward electron jet, i.e. O.lvAe- The length of the jet, 
on the other hand, becomes slightly larger for larger H . 
For Run 3, the length reached ^lOA^. 

We further examined the out-of-plane component of 
the generalized Ohm's law along x=0 as shown in the 
right panel of Figure [2] . It is expressed as 

E, = '-^^ + -VeyB^ V-r 1 

e oy c ne 

where Vg is the electron bulk velocity, T = Pexz^ + PeyzY 
is the flux of z-directed electron momentum in the re- 



connection plane (not including convection of momen- 
tum) with pe the electron pressure tensor. The observed 
electric field (solid black curve) is balanced by the sum 
(dashed black curve) of the electron inertia (red), the 
Lorentz force (brown), and the divergence of the momen- 
tum flux (green). Note that the time variation term is 
negligible and is not considered here. As is the case for 
non-moving, symmetric reconnection reported previously 
the major contribution to the inner diffusion comes 
from the pressure term. The upward outer diffusion re- 
gion is also similar to the symmetric case. On the other 
hand, the downward outer diffusion region is highly mod- 
ifled. Steeper acceleration of Vey and piled-up make 
the Lorentz term profile to be sharper. The inertia term 
with the opposite sign is also enhanced because of steeper 
gradient of Vez in the short jet. The striking feature is 
that both are enhanced in a balanced way such that 
as a whole shows the almost flat profile along the a;-axis. 
The above finding leads to the idea that the downward 
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outer diffusion region adjusts itself to buffer the effects 
of the waU such that the reconnection electric field is de- 
termined irrespective of the asymmetric boundary con- 
dition. This idea is inspected in Figure [3] which shows 
the reconnection rate, E^i, of all runs. is determined 
by averaging the electric field over a small square region 
centered at the X-line. We have verified that taking the 
time derivative of the total magnetic flux between the 
X-line and the bottom center of the pile-up region yields 
identical profile as En. For Runs 1-3 and 5, the upward 
outflow jet reaches the upper wall well after ^Icit ~160. 
For Run 4, both upward and downward jet collide at 
around ^Icit ~120. It is evident that retreating cases 
(Runs 1-3) shows almost the same profile as the non- 
retreating case (Run 4). Moreover, the profile for Run 
5 confirms the earlier result that does not depend 
on ^ [13] ■ Relatively large fluctuations are due to sec- 
ondary islands whose appearance times are indicated by 
the horizontal arrows. Also illustrated in Figure[3]are the 
timings of the jet collision with the wall as well as the 
X-line retreat. It clearly shows that the rising motion of 
the diffusion region is associated with the expansion of 
the pile- up region rather than the peak E^i- 

To summarize, the retreat motion of the X-line is con- 
stant (^O.Iva)- This speed does not depend on the ini- 



tial distance of the X-line from the wall, H, except of 
course the non-retreating case. The downward outflow 
region is largely modulated by the asymmetry. The ion 
deflection pattern depends on H and so is the length of 
the electron downward jet. The outer diffusion region 
is not elongated as in the case of symmetric reconnec- 
tion. In terms of the reconnection rate, however, it stays 
the same as the symmetric case thanks to the internal 
balance within the downward outer diffusion region. 

In addition to the initial distance from the wall, a de- 
pendence on the mass ratio /i should also be addressed. 
What we have found here is the tendency for more mag- 
netic islands to emerge with smaller fi. These islands may 
help maintain fast reconnection [ij. In fact, as shown in 
Figure [3l the fi=9 case (Run 5) generated two large is- 
lands, yet keeping the same reconnection rate as the other 
runs. Another important aspect of the /i dependence is 
the separation between the ion stagnation point and the 
X-line, but we could not measure S^s for the fj,=9 case 
because of the successive generation of the secondary is- 
lands. The mass dependence of the size of the separation 
will be studied in our future paper. Note the fact, how- 
ever, that the retreat speed was roughly equal to the ion 
flow speed so that the X-line and the stagnation will be 
collocated in the X-line rest-frame of reference. 
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FIG. 3: Dependencies of reconnection rate on the initial dis- 
tance from the wall (upper panel) and the mass ratio fi (lower 
panel). The horizontal arrows show periods of secondary is- 
land appearances. The straight vertical arrows show the ap- 
proximate times of the collision of the downward jet tcoi which 
is defined by the peak in the profiles of the downward ion fiow 
speed (see, for example, the blue arrow in Figure [T]). The 
dashed vertical arrows show the start times of the X-line re- 
treat, tret- The numbers annotated below each dashed arrow 
are the downward ion fiow maximum Vp/vA- There is no ver- 
tical arrow for Run 4 because it is the non-retreating case. 



Finally, it should be mentioned that the X-line speed 
can be faster than the self-retreating speed of O.Iva if, 
say, it is blown by another dominant X-line. Recent two- 
fluid simulations show that its reconnection rate is re- 
duced and reconnection at the moving X-line is even- 
tually terminated. Particle simulations for this kind of 
situations are also needed. 
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